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Abstract. For gravity coupled to N scalar fields, with arbitrary potential V, it is 
| shown that all flat (homogeneous and isotropic) cosmologies correspond to geodesies 

■ in an (N + l)-dimensional 'augmented' target space of Lorentzian signature (1,N), 

timelike if V > 0, null if V = and spacclikc if V < 0. Accelerating cosmologies 
correspond to timclike geodesies that lie within an 'acceleration subcone' of the 
'lightcone'. Non-flat (fc = ±1) cosmologies are shown to evolve as projections of 
geodesic motion in a space of dimension (N + 2), of signature (1, N + 1) for k = — 1 
and signature (2, N) for A; = +1. This formalism is illustrated by cosmological solutions 
, of models with an exponential potential, which are comprehensively analysed; the late- 

time behaviour for other potentials of current interest is deduced by comparison. 
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Cosmological models with scalar field matter have been much studied in the context of 
inflation and, more recently, in the context of the late-time acceleration that is indicated 
by current astronomical observations (see for a recent review). One theoretical 
motivation for these studies is that scalar fields arise naturally from the compactification 
of higher- dimensional theories, such as string or M-theory. However, the type of scalar 
field potential obtained in these compactifications is sufficiently restrictive that until 
recently it was considered to be difficult to get accelerating cosmologies in this way, 
although the existence of an accelerating phase in a hyperbolic (k = — 1) universe 
obtained by compactification had been noted [2], and non-perturbative effects in M- 
theory have since been shown to allow unstable de Sitter vacua jH]. 

In an earlier paper, we pointed out that compactification on a compact hyperbolic 
manifold with a time-dependent volume modulus naturally leads to a flat (k = 0) 
expanding universe that undergoes a transient period of accelerating expansion |E]. 
Numerous subsequent studies have shown that such cosmological solutions are typical 
to all compactifications that involve compact hyperbolic spaces or non-vanishing p- 
form field strengths (flux) El HI El 13 an d this was additionally confirmed in 
a systematic study ^Tj. Furthermore, the transient acceleration in these models is 
easily understood i n terms of the positive scalar field potential that both hyperbolic 
and flux compactifications produce in the effective, lower- dimensional, action. This 
perspective also makes clear the generic nature of transient acceleration. 

For any realistic application one would want the lower-dimensional spacetime to 
be four-dimensional, but for theoretical studies it is useful to consider a general d- 
dimensional spacetime. Assuming that we have gravity, described by a metric g mn , 
coupled to iV scalar fields (p taking values in a Riemannian target space with metric 
G a p and with potential energy V, the effective action must take the form 

J d d x V=9 (hi - l -g mn G a p{4>)d m rdn^ - V(<f>)^ (1) 

where R is the (spacetime) Ricci scalar. We are interested in solutions of the field 
equations of the action (pQ) for which the line element has the Friedmann-Lemaitre- 
Robert son- Walker (FLRW) form for a homogeneous and isotropic spacetime. In 
standard coordinates, 

ds 2 = -dt 2 + S(t) 2 dY? k , (2) 

where the function S(t) is the scale factor, and represents the (d — l)-dimensional 
spatial sections of constant curvature k. We normalise k such that it may take values 
k = 0, ±1 for a Riemann tensor Rmnop = k (gmo9np — QmpQno) with metric g on E fc . The 
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scalar fields are taken to depend only on time, which is the only choice compatible 
with the symmetries of FLRW spacetimes. The universe is expanding if dtS > 
and accelerating if d^S > 0. We need only discuss expanding cosmologies because 
contracting cosmologies are obtained by time- reversal. 

In some simple cases, the target space has a one- dimensional factor parametrised 
by a dilaton field f(4>), and the potential takes the form 



for some 'dilaton coupling' constant a and constant Vq. This model is of special interest, 
in part because of its amenability to analysis. The special case for which the dilaton 
is the only scalar field was analysed many years ago (for Vq > 0) using the observation 
that, for an appropriate choice of time variable, cosmological solutions correspond to 
trajectories in the 'phase-plane' parametrised by the first time-derivatives of the dilaton 
and the scale factor [T^J. This method (which has recently been extended to potentials 
that are the sum of two exponentials |14j ) allows a simple visualisation of all possible 
cosmological trajectories. Moreover, all trajectories for flat cosmologies can be found 
explicitly [THl HEj (see also [El EE1 HH] ) , and a related method allows a visualisation of 
their global nature [2H1 ETJ It was noted in J3| that there is both a 'critical' and a 
'hypercritical' value of the dilaton coupling constant a, at which the set of trajectories 
undergoes a qualitative change. In spacetime dimension d, these values are fHJ ITT] 



Below the critical value (a < a c ) there exists a late-time attractor universe undergoing 
accelerating expansion, whereas only transient acceleration is possible above it. The 
hypercritical coupling (a = ah) separates 'steep' exponential potentials (which arise in 
flux compactifications) and 'gentle' exponential potentials (which arise from hyperbolic 
compactification, in which case a > a c so the potential is still too steep to allow eternal 
acceleration). One aim of this paper is to generalise this type of analysis to the multi- 
scalar case. For k — 0, this has been done already for what could be called a 'multi- 
dilaton' model [22], and for the multi-scalar model with exponential potential (jSJ) [23j. 
Here we consider all cosmological trajectories (arbitrary k) for an exponential potential 
of either sign,:): and for any spacetime dimension d. In particular, we find exact solutions 
for all flat cosmologies when V < 0, following the method used in JH] for V > 0, and 
the exact phase-plane trajectories for all k when V = Vq < 0. 

| Negative exponential potentials are of relevance to pre-big bang and ekpyrotic scenarios. Their 
apparent instability can be overcome by the introduction of higher order quantum corrections to the 
gravitational action, see and references therein. 



V = V^e- 2 ^ 



(3) 




and 




(4) 
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A more ambitious aim of this paper is to determine what can be said about 
cosmologies with more general scalar potentials. What kind of model-independent 
behaviour can one expect, and how generic is the phenomenon of transient acceleration? 
Exponential potentials are simple partly because of the power-law attractor cosmologies 
that they permit, but such simple solutions do not occur for other potentials so other 
methods are needed. In this paper, we develop an alternative method of visualising 
cosmological solutions that applies to any potential, and we illustrate it by an application 
to exponential potentials. As we explain, exponential potentials serve as 'reference 
potentials' in determining the late-time behaviour of cosmologies arising in a large class 
of models with other potentials. 

Our starting point for the new formalism is the observation [23] that flat 
cosmological solutions of gravity coupled to N scalar fields with V = can be 
viewed as null geodesies in an 'augmented' target space of dimension N + 1 with a 
metric of Lorentzian signature (see for related results). Trajectories in this space 
corresponding to non-flat cosmologies are neither null nor geodesic. However, we will 
see that they are projections of null geodesies in a 'doubly-augmented' target space with 
a metric of signature (1, N + 1) for k = —1 and signature (2, N) for k = +1. 

It turns out that the extension of these results to V ^ is very simple. Flat 
cosmologies are again geodesies in the augmented target space but with respect to a 
conformally rescaled metric; the geodesic is timelike if V > and spacelike if V < 0. Null 
geodesies are unaffected by the conformal rescaling and hence continue to correspond 
to (k = 0) cosmologies for V = 0. The analysis of acceleration is also simple in this 
framework: within the lightcone is an 'acceleration cone', and if a geodesic enters its 
acceleration cone then the corresponding universe will accelerate. In the case of a flat 
target space with an exponential potential, the fixed point geodesies are just straight 
lines, and the universe is accelerating if the straight line lies within the acceleration 
cone. The further extension to k ^ cosmologies is achieved exactly as in the V — 
case. Cosmological trajectories are projections of geodesies in the doubly-augmented 
target space with respect to a conformally rescaled metric of signature (1, N + 1) for 
k = —1 and (2, AT) for k = +1. The geodesic is timelike if V > 0, null if V = and 
spacelike if V < 0. 

The plan of this paper is as follows. Section |21 reviews, in our conventions, the 
equations of motion in multi-scalar cosmologies, and we derive a useful alternative 
criterion for the acceleration of the scale factor. In section El we consider possible fixed 
points, concluding that these occur only for exponential potentials, and we determine 
the nature of these fixed points for arbitrary d, k and sign of the potential. In section 
13] we discuss the phase-plane trajectories and present some exact V < solutions. In 



Cosmology as Geodesic Motion 5 

section |S] we develop the interpretation of multi-scalar cosmologies as geodesic motion. 
In section |H1 we introduce the notion of the acceleration cone, and we show how the 
cosmologies arising in exponential-potential models fit into the new framework; other 
potentials of current interest are also considered. We summarise in section 

2. Multi-Scalar Cosmologies 

We begin by briefly reviewing the equations of motion for multi-scalar cosmologies 
that follow from the action These equations will be given in two different forms, 
distinguished by the choice of time coordinate. The second form will be useful to our 
analysis in the next section of possible 'fixed points', which correspond to power-law 
cosmologies in the first form. 

To simplify the equations we define 

(AV) a :=9t 2 f + rv/V (5) 

where is the Levi-Civita connection for the target space metric G, and we introduce 
the following 'characteristic functions' of the potential V: 

a a ((b) = ^t— 6 

2 d(j) a K J 

These are the components of a 1-form on the target space dual to a vector field a. In a 

vector notation, the scalar field equation can now be written as 

Dfo+(d-l)Hd t <f) = 2Va, (7) 

where H(t) = dtS/S is the Hubble function. The Friedmann constraint is 

S 2 [|^0| 2 + 2V-(d-l) afH 2 } = (d - 1) a~ 2 k, (8) 

where the norm \d t <p\ (and hence the inner product) is the one induced by the target 
space metric: \dt4>\ 2 = dt<fi ■ dt<fi = G a p dt4> a dt4>^ . Differentiating the Friedmann 
constraint, and using the scalar field equation (J7J), we deduce the 'acceleration equation' 

d 2 S = ^- i [a 2 c V-\d t( f>\ 2 }, (9) 

where we recall that a c is the 'critical' constant given in (0J). Clearly, acceleration is 
possible only if V > 0. Note too that H(t) cannot vanish unless either V < or k > 
(or V = 0, k = and \d t <f)\ = 0) so only in these cases can an expanding universe 
recollapse. 

For the following section, it will prove useful to rewrite the above equations in terms 
of a new time coordinate r, defined as a function of t by the relation 

dr = \V\^dt. (10) 
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Note that we have allowed for the possibility that V < 0, although V = is excluded. 
We will also set 

S(t) = e?W. (11) 

In terms of the variable (3, the condition dtS > for expansion is (3 > 0, while the 
condition d 2 S > for acceleration is 

P + (3 2 - (3 (a ■ 0) > , (12) 

where the overdot denotes differentiation with respect to r. 

In the new time coordinate, the scalar field equation becomes 

D 2 T <f> = (a- <p\ cj)- (d- l)/3(/) + 2signV r a, (13) 

and the Friedmann constraint becomes 



^\V\ 



\d>\* - (d-l)a- 2 (3 2 + 2signV =(d-l)a~ 2 k. (14) 



The acceleration equation is now 

2 • . / .\ In 2 

P = -^T 10 ' 2 " ^ + p ( a 'V + J^l sign v ' (15) 

and the condition (|T2"J) for acceleration is therefore equivalent to 

|0| 2 < assign V. (16) 

Note that this inequality is independent of k, and hence valid for all k. As expected, it 
can be satisfied only if V > 0. 

3. Fixed Points 

We now wish to determine whether the system of equations ()131 ITB]) subject to the 
constraint (fH^) . admits any fixed point solutions for which 

(3 = and D*<f> = . (17) 

We shall see that these conditions are consistent only for exponential potentials, for 
which the equations (fHfl fT3j) become those of an autonomous dynamical system, and we 
determine the fixed points and their type for all d, k and for either sign of the potential. 
For this analysis, it is convenient to introduce the quantity 



x = \i , a = |a| . (tt 



Recall that ah is the 'hypercritical' constant given in (0J). 
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3.1. Fixed Point Conditions 

Given D 2 (f> = 0, equation (JTHJ) yields 

= ca (W |0| 2 = c 2 a 2 , 3L-(p = ca 2 ), (19) 



where c solves the quadratic equation 

a 2 c 2 - (d - 1) c$ + 2 sign V = . (20) 
Given /3 = 0, equation (|TK|) yields 

2signV = (d-2)|0| 2 + ^(d-l)(d-2) [/3 2 - {3 (a • </>) 

= (d - 2) cV + ~ (d - 1) (d - 2) /5 2 - ca 2 /3 (21) 

where we have used (II 9|) to get to the second line. Using this equation to eliminate the 
2 sign V term from ()20j) we arrive at the quadratic equation 

[p - a 2 c) ({3 - a 2 c) = . (22) 

There are therefore two types of fixed point: 

• (3 = a 2 c. Using this in (|2T)|) we deduce that a ^ at and that 

e 2 = 7^. (23) 

This requires 

EITHER V > and a < a h , OR V < and a > a h , (24) 

and 

cu(j) = ±a, uj(3 = ±a 2 c . (25) 

It then follows that 

|0| 2 -(d-l)a c - 2 /3 2 + 2sign\/ = O. (26) 

This is just the Friedmann constraint for k = 0, so this type of fixed point occurs 
only for k = 0. 

• (3 = a 2 c. Using this in (|2TJ|) we deduce that 



a 



2 



c 2 = ^signV. (27) 



a 



It follows that this type of fixed point can occur only for V > 0, and that 



a, 



<£ = ±(^Ja, (3 = ±a c a. (28) 

Using this in the Friedmann constraint we deduce that 

a- 2 k = e 2 P\V\ (a 2 c -a 2 ). (29) 

We may exclude the possibility that a = a c because this yields a k = fixed point 
of the type already considered, so k ^ and we need k > for a < a c and k < for 
a > a c . Note that \<j>\ = a c for this type of fixed point, so the fixed point cosmology 
is neither accelerating nor decelerating. 
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3.2. Determination of the fixed point type 

The above results are in complete analogy with those of f3J [Ej for the one-scalar case. 
To see why, it should first be appreciated that the fixed point conditions have been 
derived on the assumption that j3 is constant, and that <p is covariantly constant. It 
then follows from ()19l l2"Uj) that a must be covariantly constant too. However, for any 
space that admits a non-zero covariantly constant vector, there exist coordinates in 
which this vector is constant, not just covariantly constant. For coordinates in which a 
is constant the potential takes the form 

V = V e~ 2a -* 1 (30) 

which is of the form (jSJ) with atp = a • <fi. Moreover, for any target space on which this 
potential is globally defined, one can find new coordinates such that 

G a/3 d<f> a dc^ = d v 2 + G (x) dx«dx? , (31) 

where G is a metric on a 'reduced' target space with (N — 1) coordinates x- m these 
target space coordinates, the scalar field equations are 

D 2 T X = a^X -(d-l)/3x, 

(p = a (p 2 - (d-l)fi<p + 2asignV , (32) 

and the Friedmann constraint is 

«c [Ixl 2 + 2 + 2 sign V ] - (d - 1) (3 2 = k (d - 1) e^'l^r 1 • (33) 

These imply the acceleration equation for (3. If we suppose that the reduced target space 
is flat then the scalar field equations, together with the acceleration equation, define an 
autonomous dynamical system: 

(p = atp 2 — (d — 1) $cp + 2a sign V , 

P= ~ T^TjT ^ 2 +lx| 2 -^signV] -/5 2 + a/30, 

X = a<px-{d-l)Px. (34) 

Note that we may consistently set x = to recover the equations of the one-scalar 
model. As there are no fixed points with ^ ^ 0, the fixed points of the multi-scalar 
model are the same as those of the one-scalar model. 

If the reduced target space is not flat then a connection term must be added to 
the last equation in (JH2J. As this introduces dependence on Xi the equations no longer 
define an autonomous dynamical system in just N + 1 variables. It is not clear how 
much difference this makes to the results: note that the connection term has no effect 
on the stability of fixed points with x — because it is quadratic in x- However, we will 
suppose for the analysis to follow that the reduced target space is flat. We may then 
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assume with no essential loss of generality that it is also one-dimensional (i.e. N = 2) so 
we have an autonomous dynamical system for just three variables. We linearise about 
a fixed point with ip — ip and $ = $ by writing 

(p = (po + x, /3 = f3 + y, X = z. (35) 

The linearised equations for (x, y, z) = U T take the form U = MU, where M is the 3x3 
matrix 

2a<p - (d - 1) $o -(d-l)<p \ 

a/5 o -4(rf-l)- 1 o a0 o -2/3 o . (36) 

a0 o -(d- 1) fa J 

The eigenvalues of this matrix determine the nature of the fixed point. We consider the 
two types of fixed point in turn: 

• k — 0. In this case 

uj(po = ±a , tufa = ±a 2 c . (37) 

The eigenvalues of M are 

± 2UJ- 1 (a 2 - a 2 ) , ±uj- 1 (a 2 - a 2 h ) , ±uj~ 1 {a 2 - a 2 h ) . (38) 

For an expanding universe we must take the top sign (corresponding, for V > 0, 
to the upper branch of the k — hyperboloid) . Then we have a stable node (all 
eigenvalues real and negative) for a < a c . For a c < a < ah we have a saddle, 
with one real-negative and two real-positive eigenvalues; the instability due to the 
positive eigenvalues is what leads to the 'flatness problem' of standard big-bang 
cosmology. 

For V < we may have a > ah, in which case the fixed point is an unstable node 
(all eigenvalues real and positive). All k — — 1 trajectories start at this unstable 
fixed point and end at the stable fixed point on the other branch of the k — 
hyperboloid; as (3 < at this other fixed point, all k — — 1 universes recollapse to a 
big crunch singularity. However, for a ^ ah there is a long quasi-static period on 
some of these trajectories. 

• k 7^ 0. In this case 

(p = ±a c , $ = ±a c a . (39) 
The eigenvalues of M are 

(d - 2) a + ay/ 16 - (d - 2) (10 - d) a 2 ] , T2a/a c , (40) 



2 



where a = 1 for one eigenvalue and —1 for another. Take the top sign again. For 
a < a c all eigenvalues are real but one is positive and the others negative, so we 
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have a saddle point. For a > a c all eigenvalues are real and negative provided either 



Otherwise (i.e. d < 9 and a > a) we have one real negative eigenvalue and two 
complex eigenvalues with negative real parts, and hence a stable 'spiral-node' (a 
spiral in the xy plane and a node in the z direction). Note that 



where the equality occurs for d = 9. Thus, for any d the k ^ fixed point is an 
unstable saddle if a < a c , and is stable for a > a c but can be either a node or a 
spiral-node. For d > 10 it is always a node, but for d = 9 it becomes a spiral-node 
at a = ah and for d < 9 it becomes a spiral- node at a = a < ah- 

In the supergravity context we are restricted to d < 11, but there is no scalar potential 
possible in d = 11 and the only possible potential in d = 10 is a positive exponential 
(with a > ah so there is no fixed-point solution). 

4. Phase-Space Trajectories 

We now have all the information needed to determine the qualitative behaviour of all 
phase-space trajectories. We begin with the one-scalar model, obtained by setting x = 0. 
The trajectories for V > are sketched in ^3] but, for the reader's convenience, we give 
a brief summary in words. All k = trajectories begin, after a big bang singularity, with 
a period of kinetic energy domination with effectively vanishing potential. For a > ah 
they end up approaching a similar late-time phase, but for a < ah they approach the 
k = fixed-point solution, which is accelerating if a < a c and decelerating if a > a c . For 
a < a c , all k ^ trajectories sufficiently close to a k = trajectory also approach the 
k = fixed point; in fact, all k — — 1 trajectories have this property. The behaviour of a 
k = +1 trajectory is determined by its relation to the separatrix of the unstable k = +1 
fixed point (which is the Einstein Static Universe for a = 0); those on one side approach 
the k = fixed point while those on the other side represent universes that recollapse 
to a big crunch. For a > a c all k = +1 universes recollapse to a big crunch whereas 
all k — — 1 universes now approach a zero-acceleration Milne universe (the k = fixed 
point now being a saddle point). For a > a the k = — 1 fixed point is an attractor spiral 
so that all k = — 1 trajectories spiral around a point of zero-acceleration; as observed in 
[Tlj . this leads to eternal oscillation between acceleration and deceleration. 

For V < the k = trajectories were studied in [21] but there seems to be no 
complete study of all trajectories for this case. We will therefore present some results 



d > 10 or (if d < 9) 



4 



(41) 



a < a = 




a c < a < ah 



(42) 
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of the one-scalar model for V < before moving on to discuss modifications that arise 
in the multi-scalar case. 



4-1. Some Exact V < Solutions 



It is well known that anti-de Sitter space is a solution of Einstein's equations with a 
negative cosmo logical constant. This space must correspond to one of the cosmo logical 
trajectories in a model with V = V < (which implies that a = 0). Set 



P = y, 



lf = X, 

to get the equations 

-(d- l)xy 



where 



x 



A = 



y = - 



y 2 + \ 2 + 



2x 



2 1 



d-1 



(43) 



(44) 



(45) 



^(d - 1) (d - 2) ' 
These equations imply that 

[x 2 - 2X~ 2 y 2 - 2] = Cx 2 '^ (46) 

for constant C. The curves with C > are the k — +1 trajectories and those with 
C < are the /c = — 1 trajectories. The two branches of the C = curve, a hyperbola, 
are the k = trajectories. The adS solution corresponds to C = — oo, and hence x = 0; 
for this special case we have 

y = -[y 2 + A 2 ] (47) 

which has the solution y = —A tan Ar, and hence 

e P = S cos Ar . (48) 

-1, as expected, and 



The Friedmann constraint implies that k 
d-l 



Vn = 



2S 2 



Ar 



Ct c t 



(49) 



and hence 



ds 2 = -dt 2 + S^ 2 cos 2 (a c f/|£b|) dT?_ 1 . (50) 

This is anti-de Sitter space. In contrast to de Sitter space, which can be realised (for 
V = V > 0) as a FLRW universe for any k (in particular as the k = fixed-point 
solution), anti-de Sitter space is realised as an FLRW cosmology only for k — — 1, and 
it does not correspond to a fixed point cosmology (at least in the sense of this paper). 
We now turn to the k = trajectories for arbitrary a. In this case we have to solve 

-(d-l)xy + a(x 2 -2) (51) 



x 
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subject to the constraint 

x 2 - 2\- 2 y 2 = 2 . 
The constraint is solved by setting 

x = v2 cosh £ , y = A sinh £ 
and the equation of motion becomes 



£ = — = [(a - a h ) e ? - (a + a h ) e ? ] . 
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(52) 
(53) 
(54) 



This equation is immediately integrated if a = ah- Having found x,y we integrate 
to deduce that 



-2a h p 



|r|e^ r 



So 



I I ~a 2 r 2 

re h 



2(d-l) 



(55) 



for integration constant So (we have set to zero the other integration constant as it can 
be absorbed into V ). 

If a ah then we introduce the new variable 



c 



ah — a 



e^ . 



(56) 



ah + a 

The two branches of the k = hyperbola correspond to positive and negative £; we take 
( > 0. We now consider separately the cases a > ah and a < a^. 

• a < ah- In this case the equation for ( is 

C = -w(C 2 + l) K C = cotwr], (57) 
where we recall that u is the function of a given in (jl8j) . Integration of (}4*3*j) now 



leads to 



'LOT COS a ah U)T. 



S [, 



sin 



a h +a 



UJT COS 



«7, 



a UJT \ 2a h" 



(5? 



a > c^. In this case the equation for £ is 

C = w (C 2 - 1) • (59) 

There is a fixed point at £ = 1, so we must distinguish between ( > 1 and ( < 1. 
One finds that 

— cothcur ( > 1 

— tanhcjr ( < 1 
For C > 1 integration of ()43|) yields 



c 



(60) 



S [smh a+ah uT cosh ah ~ a lot] ^ 



(61) 
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For ( < 1 integration of (|43j) yields 



cur cosh' 



5 [wah ah - a ur cosh a+ah lot 




(62) 



Multi-scalar Trajectories 

In the multi-scalar case, the 'spectator' fields x have an effect on the scale factor due to 
the \x\ 2 term in the acceleration equation, and in the Friedmann constraint (}3*3*j) . Phase 
space trajectories with non-zero \ wm depend on the reduced target space metric but, 
as long as there are no new fixed points, one may expect the qualitative behaviour to 
be independent of this metric. Let us therefore continue to suppose that the target 
space (and hence the reduced target space) is flat. In this case, all k = trajectories 
lie in a hyperboloid in phase space that separates the k = — 1 and k — +1 trajectories. 
This is true for either sign of V but we will now suppose that V > 0; in this case, we 
see from (fTBj) that acceleration occurs whenever the phase-plane trajectory enters the 
region with hyperspherical boundary |0| = a c . This was observed in in the context 
of an analysis of k = trajectories, but it is also valid for non-zero k. 

The only exponential models with eternally accelerating cosmological solutions are 
those for which a fixed point lies inside or on the sphere of acceleration, as happens when 
a < a c . For all other values of a, any acceleration can be at most transient, and only a 
subset of the trajectories undergo even transient acceleration; these are the trajectories 
that pass through the sphere of acceleration. The more scalar fields there are, the 
more freedom there is to avoid the sphere of acceleration, so transient acceleration in 
multi-scalar models is less generic than it is in the single-scalar model. To make this 
statement quantitative would require an understanding, along the lines of of what 
the appropriate measure might be on the space of trajectories. 

The most significant fact about a k = fixed point is whether it lies inside or 
outside the sphere of acceleration. Figure d illustrates a few examples in a model with 
two scalar fields and a positive exponential, scalar potential. 

5. Augmented Target Space and Geodesic Motion 



We now develop a method to describe the generic evolution of multi-scalar cosmologies. 
It will be convenient to define a new independent variable 7 by 

(d-l)/3 = a fc7 , (63) 
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Figure 1. These diagrams for k — cosmologies with two scalar fields moving in 
a positive, purely exponential, potential illustrate the trajectories in phase space. The 
motion takes place in the hyperboloid above the plotted plane, defined by equation $26)) . 
The constant characteristic functions a of the potential satisfy a < a c , a c < a < 
and a > ah, from left to right. The stable direction el, corresponding to j; = 0, is 
chosen at an angle 7r/4. Trajectories describe accelerating cosmologies when they lie 
inside the sphere of acceleration marked in grey. 



although we will still use j3 where this is more convenient. The Friedmann constraint 
can now be written as 

(da) 2 - \d t cf>\ 2 + ka~ 2 (d - 1) e~ 2p = 2V (64) 
and the scalar field equation as 

AV + a h {d a ) d t 4> - 2V* = . (65) 
These two equations imply the acceleration equation 

d 2 tl + — |«9,0| 2 + [(da) 2 - 2V] = . (66) 

ah y/d-l 

Using the Friedmann constraint to eliminate 2V from the scalar and acceleration 
equations we have 

AV + a h (da) dt<f> + a [|^| 2 - (daf] = k a' 2 (d - 1) e" 2/3 a , 

<9 2 7 + a h |<9t(/>| 2 = k a~ 2 a h e~ 213 . (67) 

Before considering the general case (arbitrary k and arbitrary V), we first discuss the 
special cases of V = (but arbitrary k) and k = (but arbitrary V). 

5.1. V = Case 

We may consider the (N + 1) variables 

<^ = (7,<n (68) 
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to be maps from the cosmological trajectory to an 'augmented' target space. In this 
notation, the Friedmann constraint for V = can be written as 

Gp, d^dt®" = k a' 2 (d - 1) e~ 2/3 (69) 

where 

G^ d<$>»d<$> u = -d-f 2 + G aP d<fd^ (70) 

is a Lor entz- signature metric on the augmented target space. Using the Friedmann 
constraint to eliminate the \d t <f>\ 2 term in the second equation in (JBTj) yields 

d 2 tl + a h (d tl ) 2 = -k a~ 2 (d - 2) a h e~ 2f3 . (71) 

If this is taken together with the scalar field equation then the two combine into the 
single (albeit coordinate-dependent) equation 

D 2 ^ + a h (d tl ) 8 t ^ = k a~ 2 {d - 2) a^GTd^ . (72) 

For k = this is the equation for a geodesic in a non-affine parametrisation. In terms 
of the the new time variable t' defined by 

dt! = e~ a ^dt (73) 

we have 

D 2 ^ = k a' 2 (d - 2) a h e 2{d ~ 2)p G^d ul . (74) 

For k = this is the equation of an aflinely parametrised geodesic. 

Although the cosmological trajectory in the augmented target space is not a 
geodesic for k ^ 0, it can be viewed as the projection of a geodesic in a 'doubly- 
augmented' target space of dimension (N + 2) that is foliated by hypersurfaces isometric 
to the augmented target space. Let 

tf A = (75) 

be the (N + 2) coordinates. We take the metric to be 

G AB d^ A d^ B = G^ d^d<$> v + G m (d^*) 2 (76) 

where 

= -fcajV 2 ^ . (77) 

As we now have one more variable we also need another equation. We take this to be 
the 'projection equation' 

d t ,ty* = (d- 1)2 e W-W. (78) 

With this choice, the Friedmann constraint becomes 

G AB d t ^ A d t ^ B = 0, (79) 
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and the combined scalar field and acceleration equations (|74|) are equivalent to the single 
geodesic equation 

D 2 M = 0. (80) 

Thus, cosmological trajectories for V = are null geodesies in the doubly-augmented 
target space. Note that the signature of this space is Lorentzian if k < but non- 
Lorentzian (with two 'times') if k > 0. Because \P* is not constant, the motion is not 
restricted to a single hypersurface of constant ^f* and this accounts for the fact that 
the projection of the motion onto one such hypersurface is not geodesic. This argument 
fails for k = 0, for which the projected motion is geodesic, because the metric on the 
doubly-augmented target space is degenerate when k = 0. 

5.2. k = Case 

The Friedmann constraint for k = 0, but arbitrary V, can be written as 

G^dt^dt^ = -2V . (81) 

Thus, the cosmological trajectory in this space is timelike for V > 0, null for V = and 
spacelike for V < 0. As we have seen, it is a null geodesic when V = 0. When V ^ 
the cosmological trajectories are no longer geodesic with respect to the metric G^ u , but 
they are geodesies with respect to the conformally-rescaled metric 

Gfj, u = £l 2 G^ v (82) 
v \V\e a ^ V^0 

n = < \ 2 ( 83 ) 

e a ^ V = 0. 

In fact, one finds that the equations (jBTf) are equivalent, for k = 0, to the equation 

D 2 ^ = d 2 & 1 + f» vp d t <$> u d t <$? p = fdt&* , (84) 

where T^ vp is the Levi-Civita connection for G^, and§ 

f = a h da-2a-d t <t>. (85) 

This is the equation for a geodesic in a non-affine parametrisation; noting that 

/ = dt&* dp In (2Q 2 e- ah ~<) (86) 

we deduce that the geodesies are affinely parametrised by a new time-coordinate t for 
which 

di = (2ftV Qh7 ) dt . (87) 
S It should be understood here that a = when V = as In V is undefined in this case. 



where the conformal factor is 

1 
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In other words, 

5f ^ = . (88) 

The Friedmann constraint now takes the form 

f - sign V V 4 
G,A^d^=l (89) 

Notice that t differs from t! as defined (for V = 0) by (|73|) . In fact, 

dt = 2Q 2 dt' . (90) 

This difference occurs because the afline parameter of a null geodesic is affected by 
a conformal rescaling of the metric, and the metric we are now considering is the 
conformalfy rescaled one G; clearly a null curve that is geodesic with respect to G is also 
geodesic with respect to any other conformalfy equivalent metric, such as G, but the 
affine parametrisation will differ in general. In contrast, when V 4 the cosmological 
trajectories in the augmented target space are geodesies with respect to a unique metric 
(up to homothety) in the class of metrics that are conformalfy equivalent to G, and this 
metric is G. 

5.3. General Case: Arbitrary k andV . 

The steps that led to (fHHjl for k = lead, for general k, to the equation 

{ -hafid-^e-^lVl-^dJlne^lVl) V40 
D 2 ^ = \ \ (91) 
y--ka- 2 (d-l)e- 2p G^d v {\Yie 2 P) V = 0. 

Although this is not the equation of a geodesic when k ^ 0, it is the projection of a 
geodesic in the doubly-augmented target space with respect to the conformalfy rescaled 
metric 

G AB = tt 2 G AB (92) 

where the conformal factor Q is as given in (JHHj) . and Gab is as given in (fTtjj) . From (J7HJ) 
and (|9T)j) we deduce that 

f (d - iy e~ 2f3 \V\- 1 V^O 
W = \ \ J (93) 

{(d-l)2e~ 213 V = 0. 

Given this, one then finds that (|9"T|) is equivalent to the geodesic equation 

D 2 m A = (94) 

and that the Friedmann constraint is 

, -sign\/ V 4- 
G AB d^ A d^ B ={ S r 95 

^0 V = 0. 
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Thus, the general cosmological trajectory is a geodesic in the doubly-augmented target 
space, one that is timelike for V > 0, null for V = and spacelike for V < 0. 

6. Applications 

We now consider some applications of the formalism just developed. As we are 
particularly interested in accelerating cosmologies we first consider the implications of 
the acceleration condition. We will then see how to interpret the fixed point cosmologies 
that arise for exponential potentials. This turns out to be useful when considering the 
asymptotic behaviour of more general potentials. 

6.1. The Acceleration Cone 

We first show how the condition for acceleration (JT5Jl acquires a geometrical meaning in 
the new framework. We will now assume that V > 0, as acceleration can not otherwise 
occur, and before proceeding we note that the Friedmann constraint (fbi|) can be written 

as 

7 2 = |0| 2 + 2- ka- 2 (d- \)V- X e~ 2!i . (96) 



We now introduce the following 'acceleration' metric on the augmented target space, 

1 

d~l 

and the corresponding conformally rescaled metric 



G™ d^d§ v = -- — - dj 2 + G aP d<p a d<p p , (97) 



= fi 2 G^ c , (98) 
where the conformal factor is the same as before. Noting that 

di=V2Qdr (99) 

one finds that 

grace QjpQjy = Igacc =10)2 ^ . ( 1Q 0) 

ixu t t 2 y. v 2 1 2{d- 1) 

Using the Friedmann constraint in the form (J96)) to eliminate the j 2 term, we deduce 
that 



1 



k 



GZ c di^d^ = ^ \<t>\ 2 -a 2 c +__ e -^. (101) 



2aW 



Recalling that acceleration occurs when |0| 2 < a 2 , we see that the condition for 
acceleration is equivalent, for k = 0, to 

Q»ce Q^HQ^v < . (102) 
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Geometrically, this states that a universe corresponding to a timelike geodesic trajectory- 
is accelerating when its tangent vector lies within a subcone of the lightcone defined by 
the acceleration metric on the augmented target space. 

As one might expect, an analogous result holds for k 7^ in terms of the doubly- 
augmented target space. In this case, the acceleration metric is 

GT B d^ A d$> B = G^d^d^ + — - G„ (dm*) 2 
iw {d — 1) 

= G aP d<f) a d^ + — ^— [-d^ - kafe- 2{d - 2) P {d®*) 2 } . (103) 

I & 1 ) 



(d 

Using the projection equation (JHB'j) one finds that 

|0| 2 -a c 2 l, (104) 



G a ^d^ A d^ B = \ 



where G^b = ^ 2< ^ab * s the conformally rescaled acceleration metric. The condition for 
acceleration is therefore equivalent to 

GT B di^ A d^ B <0. (105) 

This states that a universe is accelerating when the tangent to its geodesic trajectory in 
the doubly-augmented target space lies with an acceleration subcone of the lightcone. 

6.2. Trajectories for Exponential Potentials 

We now discuss how the cosmologies arising for an exponential potential of the form 
(130 j) fit into the new geometrical framework. For reasons given earlier we assume in this 
case that the target space is flat, in which case the augmented target space metric G is 
also flat, and hence G is conformally flat, with conformal factor 

Vt oc e^-^ [=► d^n oc (a h , -a)] . (106) 

Our first task is to determine the trajectories in the augmented target space that 
are associated with the fixed point solutions. Consider first the k = fixed point. This 
solution has the property that the linear function 

J($) = a 2 7 - a h a ■ (107) 

is time- independent: = J = $ M <9 M J. In other words, $ is a vector tangent to the 
hyperplanes of constant J. For N = 1 there is only one such direction, which is 
determined by the gradient of Q. Thus, $ M oc G^ v d u Vt oc {pch, a) for N = 1. In fact, this 
remains true for arbitrary iV as a direct calculation shows: 

= ±uj~ x (a h , a a ) . (108) 

A fixed point solution for k 7^ has the property that the linear function 

K{$) = a 2 c -f - a h a • <f> (109) 
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is time- independent. As surfaces of constant K are hyperplanes, a k 7^ fixed point 
trajectory is also a straight line, as a direct computation confirms: 



Thus, fixed point solutions are (particular) straight-line trajectories in the augmented 
target space. In addition, as follows from our earlier results, the k = straight-line 
trajectories are geodesies with respect to the metric G = Q 2 G, whereas this is not true 
for k ^ 0. This difference can be understood as follows: as we have assumed a flat 
target space, the metric G is flat and all straight lines are geodesies with respect to it. 
However, the relevant metric is G, and the only geodesic of G that is also a geodesic of 
G is the line of steepest descent of the function Q; as we have seen, this is precisely the 
direction of the k = fixed point trajectory. 

Note that the line determined by (jll(J|) is a generator of the acceleration cone (as 
expected from the zero-acceleration property of the k ^ fixed point). More generally, 
if a straight line trajectory lies within the acceleration cone then it corresponds to an 
eternally accelerating universe, and if it lies outside the acceleration cone (a > a c ) then 
it corresponds to an eternally decelerating universe. In the latter case it might also lie 
outside the 'lightcone' (a > a^), in which case the geodesic corresponds to the k = 
fixed point for V < 0. Note that a change in the value of the constant vector a effects 
a Lorentz transformation of the k = straight-line trajectory in the augmented target 
space. Such a transformation can take any timelike line into any other timelike line, 
and any spacelike line into any other spacelike line, but it cannot take a spacelike line 
to a timelike one or vice versa (as expected from the fact that this requires a change of 
sign of the potential). Also, it cannot take a timelike or spacelike line into a null line, a 
fact that is consistent with the absence of a fixed-point solution for a = ah- 

Generic geodesies, corresponding to generic k = cosmologies, are not straight lines 
in the extended target space, but they still have a simple description. The regime where 
the solutions are dominated by kinetic energy is given by geodesies with tangent vectors 
having large \d$</>\. This means that the generic geodesies start out in null directions, as 
could be anticipated from the fact that k = cosmologies are null geodesies when V = 0. 
Let us now follow the subsequent evolution for the V > case. As the potential becomes 
more important, the geodesies bend into the timelike cone, ultimately approaching the 
timelike straight line corresponding to the fixed point solution if a < a^. On the 
other hand, if a > ah, then the geodesic ultimately approaches a null straight-line 
geodesic. For the single-scalar case, this behaviour is shown in figure |2j Cosmologies 
accelerate precisely when their corresponding geodesies bend their tangent vectors into 
the acceleration cone. 




(110) 
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Figure 2. Generic evolution of the single-scalar, k — 0, cosmologies as geodesies on 
the augmented target space. The plot shows the {7, (p}-plane in the case of a positive 
exponential potential V with constant characteristic functions satisfying a c < a < a^. 
The decelerating fixed point solution is shown as the dashed straight line outside the 
acceleration cone. Both generic geodesies start close to a null direction when the kinetic 
energy dominates and are attracted to the fixed point direction at late times. The 
geodesic on the left goes through a regime of transient acceleration when its tangent 
vector is steep enough to lie within the acceleration cone. 



For k 7^ 0, cosmological trajectories are no longer geodesies in the augmented target 
space but they are projections of geodesies in the doubly-augmented target space. The 
vector tangent to the geodesic corresponding to a k 7^ fixed point is 

dfS A = ±-j=^ [yid^l) a, ^ a, ±y/(d=T) Vpe*^ d -W) (111) 

where we have used (193(1 to get the last entry. Note that this is not a constant (N + 2)- 
vector (even allowing for the fact that cap — (3 is constant at the fixed point) so the fixed 
point solution is not a straight line in the doubly-augmented target space, as was to be 
expected because its metric is neither flat nor conformally flat. However, using (jlllj) . 
and then (|29(1 . one finds that 

GT B d^ A d^ B = \[al~a 2 -ka; 2 \V\- 1 e-^]=0. (112) 
This again confirms the zero acceleration of the k 7^ fixed point cosmologies. 

6.3. Asymptotic Behaviour for Generic Potentials 

The asymptotic (late-time) behaviour of a large class of cosmologies in models with 
rather general potentials can be determined by comparison with the asymptotic 
behaviour of cosmologies in models with exponential potentials. To demonstrate this, 
we discuss the case of flat cosmologies with positive potentials. In this case, the scalar 
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fields must approach — > (fi at late times, where some components of the constant 
iV-vector (fio may be infinite. The characteristic functions of the potential in this limit 
are also constant: 

a = lim a(0) . (113) 

The absolute value a = |a | determines the late-time behaviour, by comparison with 
the critical and hypercritical exponents: if < a c , then there will be late-time 
acceleration, with the trajectory approaching an approximate accelerating attractor; 
obvious examples are models for which V has a strictly positive lower bound, so most 
models of interest will be those for which V either tends to zero or has a minimum at 
zero. For a c < a < ah, there is late-time deceleration corresponding to the existence 
of an approximate late-time decelerating attractor. The asymptotic geodesies in both 
these cases are timelike. In contrast, for > ah the timelike geodesic asymptotes a 
null (and hence decelerating) trajectory. 

To illustrate the above observations, we consider the following single-scalar 
examples. 

• Inverse power law potentials: 

V = l/if n . (114) 

These potentials have been much studied in quintessence models (see e.g. [Tj). 
They also arise from dynamical breaking of supersymmetry [2*%j . With <^ = oo, 
one finds a = 0. This implies late-time acceleration, even though the potential 
tends to zero at late-times. 

• That any other value of may occur is illustrated by the potential 

V = v 2m e - 2a ^ , (115) 

for which a$ = a if we suppose that if — > oo at late times. For such cosmologies 
the asymptotic behaviour is the same as it would be in the exponential potential 
model with dilaton coupling constant a. 

• It may happen that at late times the scalar fields become trapped near a minimum 
of the potential. Near such a minimum, which we may assume to occur at the 
origin of field space, the potential takes the form 

V = ip 2n + A (116) 
for integer n > 0. For positive A one finds 

a = — 117 

ip Zn + A 

which tends to zero (i.e. a$ = 0) as (p — > for any A > 0. This implies late-time 
acceleration, due to the effective cosmological constant A. However, we may now 
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consider the limit A — > 0, in which case a = oo > a^, implying that the late-time 
trajectory is null, and hence decelerating. 

7. Discussion 

In this paper we have developed a geometrical method for the classification and 
visualisation of homogeneous isotropic cosmologies in multi-scalar models with an 
arbitrary scalar potential V. The method involves two steps. In the first step, the target 
space parametrised by the N scalar fields is augmented to a larger (N + l)-dimensional 
space of Lorentzian signature, the logarithm of the scale factor playing the role of time. 
This is reminiscent of the role of the scale factor in mini-superspace models of quantum 
cosmology, especially if one views the scalar fields of our model as moduli-fields for extra 
dimensions. In this case, the Lorentzian metric on the augmented target space is the one 
induced from the Wheeler-DeWitt metric on the space of higher-dimensional metrics. 
This was the point of view adopted in [25J, where it was also observed that when V = 
all flat (k = 0) cosmologies are null geodesies in this metric (an observation that goes 
back to work of DeWitt on Kasner metrics |29j). When 7^0, flat cosmologies again 
correspond to trajectories in the augmented target space but they are neither null nor 
geodesic. However, as we have shown, all flat cosmologies are geodesies with respect 
to a conformally rescaled metric on the augmented target space, the conformal factor 
depending both on the scale factor and the potential. The conformal rescaling has no 
effect on null geodesies, of course, so flat cosmologies are still null geodesies whenever 
V = 0, but now the trajectories of flat cosmologies are timelike geodesies when V > 
and spacelike geodesies when V < 0. This is true not only for a potential of fixed sign 
but also for one that changes sign; in this case the tangent to the geodesic becomes null 
just as the scalar fields take values for which V = 0. 

The second step, which is needed for non-flat (k = ±1) cosmologies, is to further 
enlarge the target space to a 'doubly- augmented' target space foliated by copies of the 
augmented target space. One can choose the metric on this space, of signature (1, N+l) 
for k = —1 and (2, N) for k = 1 and degenerate for k = 0, such that geodesies yield (on 
projection) cosmological trajectories for any k when projected onto a given hypersurface, 
and the geodesies are again timelike, null or spacelike according to whether V is positive, 
zero or negative. This general construction thus includes all the previous ones as special 
cases. 

Accelerating cosmologies have a simple interpretation in this new framework. 
Within the lightcone (defined by the metric with respect to which the trajectories 
are geodesies) there is an acceleration subcone. A given trajectory corresponds to an 
accelerating universe whenever its tangent vector lies within the acceleration cone. In 
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particular, a flat cosmology undergoes acceleration whenever its geodesic trajectory- 
bends into the acceleration cone within the lightcone of the augmented target space. 
This can happen only if the trajectory is timelike, which it will be if the potential is 
positive. 

We have also presented in this paper a complete treatment (complementing many 
previous studies) of cosmological trajectories, for any k and any spacetime dimension d, 
for the special case of simple exponential potentials. These potentials are characterised 
by a sign (the potential may be positive or negative) and a (dilaton) coupling constant 
a (the magnitude of an N- vector coupling constant a). As has long been appreciated 
[T5] . cosmological solutions in such models correspond, for an appropriate choice of 
time parameter, to trajectories of an autonomous dynamical system in a 'phase space' 
parametrised by the time-derivatives of the scale factor and scalar fields. The qualitative 
features of these trajectories (which should not be confused with trajectories in the 
'augmented target space' parametrized by the fields themselves) are determined by the 
position and nature of any fixed points, which are of two types. There is always a k = 
fixed point unless a = ah (the 'hypercritical' value of a defined by this property) but it 
occurs only for V > if a > ah and only for V < if a > a^. There is also a k ^ 
fixed point if V > 0, which coincides with the k = fixed point when a = a c , where 
a c is the 'critical' value of a (defined as the value at which fixed point cosmologies have 
zero acceleration). 

These fixed point cosmologies have a very simple interpretation as trajectories in 
the augmented target space: they are straight lines. In the k = case these lines 
are also geodesies. For V < these geodesic lines lie outside the lightcone. For 
V > they lie inside the lightcone, but may lie either inside or outside the acceleration 
cone. For models obtained by (classical) compactification from a higher dimensional 
theory without a scalar field potential, the k = geodesic line always lies outside the 
acceleration cone, so only transient acceleration is possible in these models. For a single 
scalar field the transient acceleration is generic in the sense that it is a feature of all 
k = trajectories when a > ah and of some when a < ah- In contrast, in the multi- 
scalar case, there are trajectories that correspond to eternally decelerating universes 
even when a > ah, so transient acceleration is less generic for more than one scalar 
field. 

Although exponential potentials are of limited phenomenological value, they 
are important as 'reference potentials' in determining the asymptotic behaviour 
of cosmologies in models with other potentials, such as inverse power potentials. 
Essentially, any potential that falls to zero slower than the critical exponential potential 
(as do inverse power potentials) will lead to late-time eternal acceleration. Any potential 
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that falls to zero faster than the hypercritical exponential potential will have a late-time 
behaviour that is approximately that of a model with zero potential. 

Of course, there is no evidence for the existence of cosmological scalar fields, but 
there is evidence that the expansion of the universe is accelerating and hence for dark 
energy Whatever produces this energy, it seems reasonable to suppose that it can be 
modeled by scalar fields. If so, it is possible that future observations may be interpreted 
as telling us something about the potential energy of these fields. In view of our current 
complete ignorance of what this potential might be, we have tried, as much as possible, 
to understand generic properties, and we hope that the methods developed here will be 
of further use in this respect. 

Note added 

The interpretation of solutions of gravitational theories as geodesies in a suitable metric 
space has a considerable history of which we were mostly unaware at the time of writing 
this paper. The geodesic interpretation of k = cosmologies that we have presented 
is closely related to the Maupertuis-Jacobi principle of classical mechanics. Relatively 
recent work on this topic includes |3*m IBTj. 132*] . A particle mechanics formulation of 
the geodesic interpretation described here, including our proposed extension to k ^ 0, 
has recently been found, and applied to a model in which cosmological singularities 
correspond to horizons in the augmented target space [33] . 
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